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Abstract 

Milgrom's Modified Newtonian dynamics (MOND) provides an efficient way to summarize phe- 
nomenology of galaxies which does not lean on the notion of dark matter; it has great predictive 
power. Here I briefly review MOND as well as its implementation as a nonrelativistic modified 
gravity theory, AQUAL. Gravitational lensing and cosmology call for a relativistic gravity theory 
different from general relativity if dark matter is to be avoided. In recent years such a theory, 
TeVeS, has emerged from the marriage of AQUAL with the timelike vector field of Sanders. I 
discuss its structure and some of its successes and shortcomings. 
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1. Introduction 

As other talks in this conference attest, the dark matter (DM) paradigm is a dominant 
one in physical science, on Earth and away from it. So why look for alternatives to 
it? There are three good reasons. First, although many astronomical observations seem 
to demand the existence of DM in large quantities, this is always inferred through the 
standard law of gravity. Thus if our understanding of gravity on astronomical scales 
is flawed, that inference is at grave risk. Second, apart from sporadic reports which 
often clash with known constraints, DM has not been identified directly, i.e., by means 
of a nongravitational experiment. Third, history suggests that it is best to arrive at 
the accepted picture of a phenomenon by confronting rather different paradigms for it; 
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confining all attention to the DM paradigm is dangerous until we understand more about 
the problem it is supposed to solve. 

In this review of the MOND alternative to DM, I proceed from some empirical facts 
about galaxies to the MOND paradigm, and thence to modified gravity theories. 

Disk galaxies rotate; the stars and gas move together in circular concentric orbits. The 
circular velocity is best deduced from the Doppler shift of radio waves coming from atomic 
hydrogen gas, or that of optical lines from the stars. When plotted against radius r, this 
velocity constitutes the rotation curve of the galaxy. According to a simple, Newtonian, 
balance of centrifugal acceleration and gravitation pull of the galaxy, the rotation curve 
should first rise (because more mass is being included with increasing radius) and then 
decrease as i/^/r . But the rotation curves of real galaxies, after first rising, almost 
invariably flatten out, and this behavior continues well beyond the edge of the optical 
disk as traced by the 21 cm line of atomic hydrogen. The reticence of rotation curves to 
fall off as theoretically expected counts as one of the principal reasons for the widespread 
belief that disk galaxies are all embedded in massive DM halos. 

And what sets the velocity in the flat part of the rotation curve? A simple rule, the 
famous and extremely useful TuUy-Fisher law. In one modern version of it, the infrared 
(K-band or 2.2 micrometer) luminosity of a disk galaxy is proportional to the fourth 
power of the circular velocity in the flat part ■ l] . In its "baryonic TuUy-Fisher" , version 
the mass in baryons (stars and in gas) is proportional to the fourth power of the circular 
velocity; specifically M/Mq « 50(u/km s^^)** [2 . In brief, the velocity is tied directly to 
the visible matter content. 

In the DM halo paradigm such a law must follow from considerations of galaxy forma- 
tion, a very messy process by all accounts. It has been hard to understand the sharpness 
of the correlation in this view. And the specific form of the power-four law is not well 
understood either. 

In the large clusters of galaxies, agglomerations of hundreds of galaxies moving with 
random velocities of up to 10"^ km s~^ through a medium of very hot gas, Newtonian virial 
theorem estimates of a cluster's mass from the galaxy velocities are very large compared 
to the mass seen directly as galaxies and gas. This mass discrepancy also shows up when 
the overall cluster mass is determined Newtonially from the assumption that the hot gas 
is in a hydrostatic state, or when gravitational lensing by a cluster is analyzed in the 
framework of general relativity (GR). The conventional solution is to assume that the 
typical cluster contains DM to the tune of about five times the visible mass [T]. 

We shall set c = 1 in most equations. 



2. The MOND paradigm and AQUAL 



Milgrom's MOND paradigm [314J grew out of dissatisfaction with the DM idea. Mil- 
grom replaces the standard relation between acceleration of a test particle and the am- 
bient Newtonian gravitational field, ma — — V"I>n, by 

/i(|a|/ao)a = -V$N- (1) 

where $n is ascribed to the visible matter alone. Here qq ~ 1.2 • 10^^°ms~^, a preferred 
acceleration, is of the order of the centripetal accelerations of stars and gas clouds in 
the outskirts of disk galaxies. The il{x) is a positive function with // — s- 1 when a 3> 
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ao, so that we get back Newtonian behavior as appropriate in the solar system where 
\d\ ~ 10~^ms~^, and fl{x) « x when a <C ag. This last is the deep MOND limit which 
should be relevant in the outer parts of disk galaxies. If for it we set |a| = v'^/r and 
|V$n| — GM/r^ {AI total galaxy mass), we find the MOND formula to be satisfied for 
V =const. Thus in the outskirts of a disc galaxy, where mass no longer grows with radius, 
MOND predicts a flat rotation curve. And from the coefficients in the above we get that 
M = v'^/Gao, which is just the baryonic Tully-Fisher law. Because the mass to infrared 
luminosity is fairly constant among disk galaxies [1] , the infrared luminosity version of the 
Tully-Fisher law follows too. One can view oq as determined by the observed Tully-Fisher 
law. 

How well does the MOND paradigm reproduce the shapes of observed disc galaxy 
rotation curves in terms of the mass distribution actually seen in them (gas and stars)? 
Fits to over one hundred disk galaxies bear witness that MOND is very successful |ll5l6j . 
With the choice fl{x) = x{l + x)~^ (the simple function [7]) and with the mass-to- 
luminosity ratio T for the stellar population as its only fitting parameter, MOND is 
as successful as DM halo models with three fitting parameters (typically T, central 
velocity dispersion of DM particles and scale of length of halo). In addition, the features 
in a MOND predicted rotation curve correlate well with features in the visible mass 
distribution of the relevant galaxy. This is embarrasing for the DM paradigm in which 
the gravitational field that controls the rotation curve comes mostly from the DM. It 
is no longer controversial that for disk galaxies of all scales, MOND provides the most 
economical description. 

While the MOND equation ([T]) is incompatible with the conservation laws (except if 
applied solely to test particles), the MOND paradigm can be made dynamically consistent 
by reformulating it as a modified gravity theory, AQUAL, derived from an AQU adiaXic 
iagrangian [8] . In AQUAL the Poisson equation is replaced by 

V-[/i(|V$|/ao)V$] =47rGp, (2) 

where p is the baryonic matter mass density, and — gives the acceleration field. For 
spherically symmetric systems the first integral of the AQUAL equation with —a 
reproduces Eq. ^ exactly. For disk like systems the two differ by some 10-15% [9]. 
Thus given the notoriously inaccurate character of astronomical data, AQUAL fits disk 
galaxies just as well as does the MOND equation. 

Not only is AQUAL equipped with the usual set of conservation laws, it also removes 
the MONDian paradox of why stars whose entrails are mostly strongly accelerated ions 
move in a galaxy in non- Newtonian fashion [8J. AQUAL has established [8] the long 
conjectured external field effect [3^: the suppression of MOND efl^ects in a system im- 
mersed in a sufficiently strong gravity field. And it has permitted full understanding of 
how MOND stabilizes galaxies against the "cold disk" instability [lOlllj which supplied 
the original raison d'etre for DM halos. Today AQUAL is the basis for several A^-body 
simulators which are elucidating the role of MONDian effects in galactic formation and 
evolution [12|13j . (For more on MOND and AQUAL see recent reviews, e.g. [14|15j ). 

But AQUAL is powerless to give a complete account of the dynamics of the large 
clusters of galaxies. In one of these the accelerations of the galaxies near the center can 
well exceed ap, so both MOND formula and AQUAL would predict nearly Newtonian 
behavior for cluster cores. Thus while MOND- AQUAL models of cluster dynamics explain 
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some of the mass discrepancy, they leave a factor of 2-3 unexplained [T]. Unlike galaxies, 
clusters seem to require invisible mass. The promising idea |16| that the invisible mass 
is made up of neutrinos (which are known to exist) has been running into problems 
recently [17] . Perhaps clusters contain large amounts of dark baryons [18] . Or the problem 
might be solved in MOND spirit if a natural way were found to shift upwards Milgrom 
acceleration scale oq in systems bigger than galaxies. 

Nonrelativistic by construction, AQUAL cannot serve to describe gravitational lensing 
by galaxies or clusters of galaxies, or cosmology, both essentially relativistic in nature, 
and strongly bound up with the DM issue. In fact, just a few years ago it was common for 
DM pundits to reject MOND because "it cannot be framed relativistically" . There were 
indeed various stumbling blocks on the way to this goal, but they proved surmountable. 
Today there are an handful of relativistic MOND gravity theories which disprove that 
pessimistic assessment [191201211 . We turn to the earliest and perhaps simplest of these. 

3. Tensor- Vector Scalar theory and gravitational lensing 

The requisite relativistic theory of gravity must depart from GR — the standard gravity 
theory today — because the latter has a uniformly Newtonian nonrelativistic limit. But 
it should not depart from it too strongly because GR is known to be accurate at solar 
system scale. GR embodies the equivalence principle; in its weak form this principle 
requires that all matter propagate on the same metric. Because of its philosophical basis 
and strong support from Earth-bound experiments, we would like to retain that feature. 
This means all matter actions in the new theory must be written with one metric, gaf3- 
If we write the gravitational action in Hilbert-Einstein form using ga/3 to form the scalar 
of curvature, R, we are back to GR (unless we are willing to introduce new types of 
gravitational fields). We may choose to depart from GR by modifying the gravitational 
action, e.g., as in f{R) theories. But whatever successes maybe claimed by these last 
in regard to the dark energy puzzle, they cannot reproduce the regularities shown by 
galaxies and summarized by MOND. 

By contrast, TeVeS retains the Einstein-Hilbert gravitational action, but written with 
a metric gap distinct from gaf3] this gives rise a controlled departure from GR. Relating 
gap and gap by a conformal factor is a simple possibility. This leads to a relativistic 
version of AQUAL f8]. However, the scalar field is found to be afflicted by superluminal 
propagation, and the theory does not provided for the enhanced gravitational lensing 
that would be required to dispense altogether with DM in galaxies. 

A significant idea here is Sanders' unit timelike 4- vector field Ua Instead of the 
conformal relation between the metrics, which is the root of relativistic AQUAL's failure 
in the gravitational lensing domain, Sanders used 



{(j) is a dimensionless auxiliary scalar field), which can give lensing of the correct strength. 
But since Sanders regarded Ua as constant and pointed in the time direction of the 
cosmos, the resulting theory was not covariant. 

TeVeS [19j achieves covariance by converting Ua into a dynamical field with action 



g^p = e-^^gap - {e^^ - e~^^)UaUp 



(3) 
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where K ia dimensionless parameters and A is a Lagrange multiplier field charged with 
maintaining the normalization of lAa to —1, which is the easiest way to make the vec- 
tor everywhere timelike. (Of late it has been found useful to buttress with a term 
K{g°^^Ua-,l3Y to prevent the formation of caustics of the integral lines of g°'^Up and 
inconsistencies in the weak field approximation to the theory). 
The action governing this scalar's evolution is 

Ss = -^J^ J T{ke{g^P -g^^g^'U,Us^,^<j^.,p) {-gf'd'x, (5) 

where J- is a, positive function, k is another dimensionless coupling constant, and ^ is a 
constant scale of length. This differs from a covariant version of the AQUAL lagrangian [8] 
only in the appearance of g°'^ — W^U^ instead of g°'^ , a change introduced to forestall 
superluminal tp propagation. In the limit fc ^ with K (x k and i oc /c~^/^, TeVeS 
reduces to GR jT9]. More on TeVeS can be found in the reviews [14|23j . 

For weak $n and (p and a quasistatic situation, TeVeS yields the linearized line element 
(corresponding to gafs) 

ds^ = -(1 + 2$N + '2(p)dt^ + (1 - 2$N - 2(j)){dx^ + dy^ + dz^). (6) 

The scalar (/)'s equation takes on the form of Eq. [2] sourced by the same p that sources 
<i>N, with p{y) = dT{y)/dy playing the role of fi. With the choice p,{y) = D^y for 
< J/ <C 1 (-Da positive constant), and /i(j/) — > 1 for y oo, the role of ao is played by 
c'^^/k{A7TD£)~^ Whenever |V(/)| ^ ao, the scalar dominates over $n in the metric, and 
we have AQUAL behavior of dynamics with $ w 0. And whenever |V(/)| ^ ao, 4> has the 
same form as ^n, so that TeVeS reduces to Newtonian gravity, albeit with a rescaled 
gravitational constant. In these senses TeVeS is a relativistic MOND gravity. 

According to Eq. [SI in linearized approximation in TeVeS, just as in GR, one and the 
same potential rules dynamics and controls gravitational lensing. If TeVeS describes well 
the dynamics of a system everywhere, it should predict the same pattern of gravitational 
lensing that GR would predict, were the latter supplemented with the appropriate dis- 
tribution of DM to fit the dynamics everywhere. In practice this obvious dictum may 
work less than perfectly when the extant data pertains to just a limited part of the sys- 
tem. TeVeS does well in the confrontation with data for strong lensing (multiple images) 
by galaxies |24I23I14] . It meets problems for weak lensing (distorted single images) by 
clusters of galaxies |25|26|23] : basically it does not account for the observed distortion 
without the help of invisible matter apart (in addition a reasonable dose of neutrinos). 
This problem is dramatized by a handful of colliding clusters, and has been made much 
of. But lamenting over TeVeS would be premature. MOND has never dealt perfectly 
with the dynamics of clusters. So TeVeS, which was designed with MOND in mind, 
could not expected to do well in this business, and modifications of it may be in order. 
And, as already mentioned, clusters may well contain large amounts of as yet invisible 
matter 18J. Enlightenment from the MOND paradigm is not at an end! 
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